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We consider associated production of squarks and gluinos with the lightest supersymmetric particle (LSP), 
or states nearly degenerate in mass with it. Though sub-dominant to pair production of color S[/(3) -charged 
superpartners, these processes are directly sensitive to the wavefunction composition of the lightest neutralinos. 
Exploiting event-shape variables - including some introduced here for the first time - we are able to identify the 
composition of the LSP by selecting events involving a single high-pr jet recoiling against missing transverse 
energy. We illustrate the proposed technique on a set of benchmark cases and propose methods for applying 
these results in more realistic experimental environments. 



I. INTRODUCTION: GENERAL GOALS 

With the LHC experiments currently collecting data it is 
not inconceivable that a discovery of low-energy supersym- 
metry could be made in the near future. The subsequent 
months and years will be spent making numerous measure- 
ments of masses, cross-sections and branching ratios. As was 
emphasized in earlier work for the theorist trying to recon- 
struct the underlying supersymmetric Lagrangian and looking 
for clues as to the oiigin of supersymmetry breaking, this in- 
formation is useful only to the extent that it can (uniquely) 
determine the soft supersymmetry breaking parameters them- 
selves. Of these, it has been argued that the soft masses of 
the gauginos are particularly important for distinguishing be- 
tween high-scale models of supersymmetry breaking H. In 
recent work a subset of the authors demonstrated how a syn- 
thetic approach that considers an ensemble of targeted observ- 
ables can be used to perform a fit to the ratios of gaugino 
masses Here we would like to turn our attention to the 
other great unknown of the gaugino sector: the composition 
of the wavefunction of the lightest neutralino. 

It has long been appreciated that in the R-parity conserv- 
ing minimal supersymmetric Standard Model (MSSM), the 
lightest supersymmetric particle (LSP) will be stable and can 
therefore provide a good dark matter candidate if it is un- 
colored and electrically neutral |j4l|. Assuming that this LSP 
is the lightest mass eigenstate of the neutralino system, then 
the cosmological properties of this state will be highly depen- 
dent on the composition of its wavefunction. For example, the 
question of whether the thermally-produced reUc density of 
the LSP is sufficient to account for the known non-baryonic 
dark matter - or whether some non-thermal production mech- 
anisms will be needed - is crucially sensitive to the wavefunc- 
tion of the LSP jsi]. So too is its scattering rate in terrestrial 
direct detection experiments J^. The question becomes even 
more urgent if the measured excess in the positron composi- 
tion of cosmic rays ||7[ Is]] is taken to be a signal of new (su- 
persymmetric) physics igl lToll . 

Extracting the wavefunction of the LSP from measurements 
of cross-sections x branching fractions is known to be ex- 
tremely difficult at hadron colliders such as the LHC due to 
an inability to over-constrain the independent entries in the 
neutralino mass matrix ifTTll . The typical strategy therefore in- 



volves a global fit to these entries, which typically suffers from 
a lack of uniqueness (the so called 'inverse problem') lfl2ll . 

In this paper our goals are more modest: we will seek a set 
of observables which are directly sensitive to the wavefunc- 
tion composition of the LSP without attempting to reconstruct 
the full set of eigenvalues and eigenvectors for the neutralino 
system. The analysis techniques presented here are to be un- 
derstood as a suggested first step towards a realistic measure- 
ment strategy for this important theoretical property. Here we 
will demonstrate the effectiveness of our analysis technique 
on test cases based on the highly-studied Snowmass bench- 
mark points ifTsIl . We expect that the techniques we will ex- 
plore will prove fruitful only after a sizable period of data- 
taking, and we therefore consider the case of O{l0) fb^^ of 
integrated luminosity at =14 TeV center of mass ener- 
gies. 

In Section we present a non-technical summary of the 
basic approach and introduce the benchmark models we will 
employ as illustration. In SectionHIIlwe introduce a set of new 
event-shape variables which will prove crucial to the analy- 
sis method which follows. We will use the benchmarks to 
demonstrate the correlation between LSP wave-function ex- 
tremes (bino-like versus wino-like versus Higgsino-like) and 
the event-shape distributions for associated production of a 
single LSP with SU (3)-charged superpartners. The overall 
rate gives information on the Niqq coup ling , which is sen- 
sitive to the wavefunction of the LSP. In |14{| it was demon- 
strated that by looking for supersymmetric monojets, the cou- 
pling may be determined to (10%) for a large region of wino- 
like LSP parameter space. 

In Section HV] we provide the majority of our results. We 
will demonstrate how extraction of events with direct produc- 
tion of the LSP, in association with a strongly-coupled super- 
partner, can be performed with high efficiency from the gen- 
eral superpaitner production modes. It will be shown that a 
combination of distribution shapes and integrated count rates 
in select channels can distinguish between wavefunction ex- 
tremes, when other aspects of the superpartner spectra are held 
fixed. We will indicate how the overall signal can be separated 
from the Standard Model background at a cursory level - a 
more detailed treatment of backgrounds will be reserved for 
a future analysis. Some directions for improving the analysis 
will be given in the concluding section. The analysis that we 
are discussing will be done after superpartners are discovered. 
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and their masses are known and can be used to simplify the 
studies. 



n. OVERVIEW OF THE IDEA 

In the MSSM the neutrahno sector consists of four states 
whose masses are given (at tree level) by the eigenvalues of 
the neutralino mass matrix 



/ Ml 


-swCpMz 
\ swspMz 





Ms 
cwcpMz 
-cwSpMz 



-swCpMz 
cwcpMz 





where Mi is the soft supersymmetry breaking mass of the hy- 
percharge U(l) gaugino at the electroweak scale, M2 is the 
soft supersymmetry breaking mass mass of the SU(2) gaugi- 
nos at the electroweak scale, cw = cos 6w and sw = sin 6w 
involve the weak mixing angle, and cp ~ cos/3 and sp ~ 
sin /? involve the ratio of the two Higgs scalar vevs (tan /3 ~ 
(hu) I {hd)). The matrix O is given in the (B, W, H^, Hi) 
basis, where B represents the bino, W represents the neutral 
wino and _ff ° and are the down-type and up-type Higgsi- 
nos, respectively. 

The mass matrix in ([T]i is diagonalized by a unitary matrix 
N whose entries we will denote by Nij. The content of the 
LSP can therefore be parameterized by the expression 

iVi = NiiB + NviW + A^iaffS + ^i4-ff°, (2) 



1. Thus in 



which is normalized to Nf^ + Nf^^ 
saying that the bino content of the lightest neutralino is high, 
we mean Nl^ ~ 1, for a wino-like LSP we mean — 1 
and similarly for the Higgsino limit we mean N^^ + N^j^ ~ 1. 

Ideally we would like to find a set of observables which 
will accurately measure each of the entries Nij without at- 
tempting to reconstruct the entire matrix N. In this paper we 
will concentrate on distinguishing between the three extreme 
cases defined in the previous paragraph, reserving a more gen- 
eral treatment (with arbitrary mixtures of various component 
states) to a future study. 



A. Relevant Processes 

To study the wavefunction of the LSP at the LHC it will be 
necessary to isolate processes that depend strongly on the en- 
tries of the eigenvector Nij. We therefore consider associated 
production of the lightest neutralino with other superpartners 
- particularly the squarks and gluinos, which we will refer 
to as 'semi-strong' production modes. Associated production 
with other neutralinos and charginos may also be important, 
and will be discussed in more detail below. Table |T] lists the 
allowed production processes xNi we will consider for the 
extreme cases for our neutralino (the limit of 100% bino or 
wino or Higgsino content). A checkmark indicates that this 
process is allowed for the particular wavefunction extreme. It 
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TABLE I. Allowed production processes of the form x Ni, for the 
pure wavefunction limits. A checkmark means that process is al- 
lowed. In this table q always represents a squark of the first two 
generations and we assume no mixing between the supeipartners of 
the left- and right-handed quarks. 



Channel 


a(fb) 


Ng 


88.4 


NqR 


219.3 


NqL 


18.2 


NN2 


0.9 


NCi 


2.7 



TABLE II. Production cross sections for processes of Table IDfor the 
particular example of SPS point lA. The total supersymmetric pro- 
duction IS C7suSY — 41.5 pb 



is immediately clear from the 'texture' in Table |l] that these 
processes carry sufficient information to distinguish between 
the extreme cases. Naturally as one interpolates between ex- 
treme cases the entire table is filled out to varying degrees, but 
here let us focus on the simpler case of understanding the na- 
ture of the pure LSP limits. Focusing on the first three lines 
(the semi-strong production processes), which generally have 
higher cross-sections than electroweak gaugino pair produc- 
tion, the overall strategy is apparent. If we can make mea- 
surements so as to fill in this triplet with one's and zero's we 
would match it to the patterns in Table U and identify the com- 
position of the LSP. 

There are, however, a number of issues which complicate 
matters. The processes in Table |I]typically have cross-sections 
which are at least an order of magnitude smaller than domi- 
nant processes like gluino pair production or gluino/squark 
associated production. Thus we expect that the measure- 
ments we hope to make require a significant amount of in- 
tegrated luminosity. For example, the much-studied Snow- 
mass point SPS 1 A has a mostly bino-Hke LSP (Nf^ = 0.996) 
and relatively light gluino (rrig ^ 600 GeV). It is therefore 
an excellent candidate for early discovery at the LHC. Us- 
ing PYTHIA 6 . 4 lfl9ll we calculate the total supersymmet- 
ric production cross-section to be CTsusy = 41.5 pb, almost 
half of which comes from gluino/squark associated produc- 
tion (jgq = 20.6 pb. The production cross-sections for the in- 
dividual sub-processes of Table |T] are given for SPS point lA 
in Table HIl Note that these values are given in femptobarns. 
Note that PYTHIA by default introduces no mixing between 
the superpartners of the left-handed and right-handed quarks 
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Extreme 


Number 


Identities 


Pure Bino 


1 


iVi 


Pure Wino 
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iVi,Ci 


Pure Higgsino 
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Ni,Ci,N2 



TABLE in. Effective number of LSPs producing ]^t, for the extreme 
wavefunction endpoints. Note that the Higgsino extreme may even 
have four effective LSPs if we include the N3, which is often also 
close in mass. 



Base Model 


mo 


■mi/2 


Ao 


SPS lA 


100 


250 


-100 


SPS lA' 


1000 


250 


-100 


SPS 2 


1450 


300 





SPS 2' 


200 


300 






TABLE IV. Input parameters for the four base models. All masses 
are given in units of GeV and all models have tan /3 = 10 and fi > 0. 



for the first two generations of squarks, which are here de- 
noted by the generic symbol q. Extracting these subdominant 
processes from the supersymmetric and Standard Model back- 
grounds will be challenging, but not impossible, as we will 
demonstrate in Section HV] 

An additional, and more subtle, challenge will be isolat- 
ing only those cases in which the lightest neutralino A^i is 
produced in the semi-strong production process, as opposed 
to heavier gaugino states. To illustrate, consider the case of 
an extremely wino-like LSP. To achieve this outcome it is 
typically necessary to have M2 ^ Mi, jj. in the neutralino 
mass matrix. This simultaneously produces a light chargino 
Ci which is then typically nearly degenerate in mass with 
the lightest neutralino. The phenomenology of this extreme 
case is very different from the more familiar bino-like extreme 
for which rrij^^ ~ m^^ with both being significantly more 
massive than the LSP. The wino-like extreme has been stud- 
ied extensively in the context of models of anomaly-mediated 
supersymmetry breaking ll20l - j22ll where there are effectively 
two LSPs since the decay of the lightest chargino to the light- 
est neutralino involves very soft decay products which are 
generally not detected. The same analysis can be performed 
for the case of the Higgsino extreme, for which we require 
ji ^ Ml, M2 and often have three, or even four effective 
LSPs quite close in mass with one another This state of af- 
fairs is summarized in Table |III] In the Higgsino case this is 
especially vexing as the gauginos A^2 and Ci may very well 
appear in semi-strong production processes, even when the 
true LSP A^i cannot. As we will investigate in Section lTVl this 
wiU turn out to be one of the largest difficulties in extracting 
the wavefunction of the LSP from LHC data. 



B. Benchmark Models 

To focus our study we will work with a set of bench- 
mark models, sacrificing some generality for concreteness. 
Our starting point will be the pair of benchmarks SPS lA 
and SPS 2 from the Snowmass benchmark set ifisll . These 
two "base models" are derived from minimal supergravity 
(mSUGRA) which postulates an overall scalar mass mo, over- 
all gaugino mass mi/2 and overall scalar trilinear coupling 
at some high energy scale, which then must be evolved to 
low energies via the renormalization group (RG) equations. 
The points SPS lA and SPS 2 are designed to yield oppo- 
site hierarchies between the lightest 5t/(3)-charged super- 
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300 


815 


101 
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211 
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711 




Wino 


300 


98 


815 


103 


207 




99.0% 


654 


711 




Higgsino 


387 


815 


108 


103 


13 
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98.1% 


654 


719 


SPS 2 


Bino 


98 


300 


815 


101 


214 


214 


99.6% 


783 


979 




Wino 


300 


98 
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104 


207 




99.0% 


783 


979 




Higgsino 
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815 
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104 


13 
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98.2% 


783 
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98 


300 


815 
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206 


206 


99.6% 


714 


482 




Wino 


300 


98 


815 


101 


204 




99.0% 


714 


482 




Higgsino 


400 


815 


108 


103 


13 


6 


98.1% 


715 


503 



TABLE V. Input Lagrangian masses and physical eigenstate masses 
for the twelve benchmark points we will consider in what follows. 
The values of A = mjva — mjvj and A* — mci — mjvj reflect the 
degeneracy in effective LSPs shown in Table |llll The LSP mass IS 
kept the same for all the models, so kinematical issues do not confuse 
the analysis. 



partners. Thus SPS lA has rrig > m^lniin while SPS 2 has 
iTig < mglniin, with the lightest squark being a stop in both 
cases. As we will soon see, heavier squarks will have strong 
implications for the methods we develop in the next section. 
With this in mind we will develop a variant of each of these 
two points with the opposite gluino/squark ordering. The set 
of mSUGRA input parameters for the four benchmarks are 
given in Table [TVl with /i > and taii/3 = 10 for all four 
points. 

Over a vast amount of the allowed parameter space of the 
mSUGRA paradigm the lightest supersymmetric particle is 
an overwhelmingly bino-like neutralino. The models in Ta- 
ble IIVI therefore, provide a good array of superpartner mass 
patterns but absolutely no variety in the nature of the LSP 
wavefunction. To remedy this deficiency it will be necessary 
to modify the input parameters of the neutralino sector from 
those derived from the values in Table |IV] We will do this 
in such a way as to keep as much of the particle spectrum 
fixed as is possible, particularly the mass of the lightest eigen- 
state. Specifically we will choose sets of the input parame- 
ters {Ml, M2, fJ,} which achieve at least 98% purity for each 
of the three wavefunction extremes: bino-like, wino-like and 
Higgsino-like. 

For the bino-like extremes we use values very close to those 
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that arise from the RG evolution of the mSUGRA inputs in Ta- 
ble |IV] and subsequent requirement of electroweak symmetry 
breaking: Ah = 98 GeV, Ma = 300 GeV and n = 815 GcV. 
To achieve the wino-like extreme in all four cases it is suffi- 
cient to exchange the values of Mi and A/2 at the electroweak 
scale. This roughly maintains the mass of the LSP as well 
as the mass difference A = thin^ — m between the light- 
est neutralino and the second lightest eigenstate. Note that this 
exchange has no effect on the gluino or scalar fermion masses, 
but the wino-like limit always implies a mass difference be- 
tween the lightest chargino and the LSP which is vanishingly 
small = —rriM^ < 1 GeV. To achieve the Higgsino- 
like extreme we set M2 = 815 GeV at the electroweak scale 
and then choose the values of Mi and ^ according to the ap- 
proximate formula 



^(Afi) ^ (3.3 X 10" 



M- 



-2.75\ 



-1/2.7 



(3) 



in such a way as to keep the LSP mass constant and in agree- 
ment with the bino-like and wino-like extremes. In this case 
both A and A* are required to be small but non-vanishing. 
These values are collected in Table [V] Note that changing the 
value of the //-parameter affects the squark and slepton masses 
through the off-diagonal mixing terms, but the differences be- 
tween the masses of these states for each triplet of models is 
small. 

Finally, let us note that the rather large production cross- 
sections for supersymmetry associated with point SPS lA 
may imply that this precise model point has already been ex- 
cluded by ATLAS and CMS searches for events with jets plu s 
missing transverse energy in the first 35 pb^^ of data lUSUlTll . 
This need not be true, however, of the variants with implied 
non-universalities in the gaugino sector which we have con- 
structed and listed in Table IV] [fisll . Therefore, given the thor- 
ough study these 'standard candle' models have received in 
the literature, and as we will not be discussing discovery of 
supersymmetry in this work, we will continue to work with 
the set of models given in Tables HVlandrV] 



III. EVENT SHAPE VARIABLES 

Looking at the semi-strong processes in Table I] we see 
that the production channels that will interest us are ones that 
will have unbalanced visible energy in the rest frame of the 
primary collision. If we restrict our attention, therefore, to 
the transverse plane, event variables which capture this lop- 
sided nature should be helpful in distinguishing these pro- 
cesses from the much more dominant gg and gq production. 

A widely used and very familiar quantity associated with 
these types of event-shape variables is sphericity, s, which is 
defined by lEl 



PaiPbi 



x,y,z, 



s^l^ (4) 



2 Tr(5) 



where Ai < A2 < A3 are the eigenvalues of the matrix S. 
When restricted to the transverse plane the relations in (2) be- 



nab 



PaiPbi , a. 



b = x,y , ST = 4 



Dct(SV) 
Tr(5T)2 



(5) 



For the processes that will interest us these variables have very 
similar distributions across all dominant SUSY production 
channels. It is for this reason that cuts on transverse sphericity 
are often imposed in inclusive analyses that involve multijet 
events ll24ll . 

To get at our semi-strong processes, therefore, we will need 
to look beyond the sphericity variables. The next class which 
will prove useful for our purposes are the recoil variables, r, 
which are related to observables such as thrust. In this paper 
we will utilize a triplet of such variables, defined by 
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eT^ 



(6) 



(7) 



(8) 



Ti\ 



In addition to the above, we will also introduce a set of new 
variables, q. 
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(Ei i\) i,j<i 



J2 ^ 



(9) 



qr 



2/7r 



(Eji^T.ir 



^ I^T.|I^T,||(0»-0,)|. (10) 
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We will use the ensemble of twelve models in Table [V] to 
demonstrate the efficacy of our shape variables in separat- 
ing semi-strong production processes from one another and 
from the remainder of the supersymmetric production pro- 
cesses and the Standard Model backgrounds. We begin by 
looking at the transverse objects tt, r'j, and qT for the origi- 
nal SPS lA benchmark (the first model in Table IVTi. 

The distributions of these variables are given in Figure [1] 
The shape variables are formed from the momenta of all vis- 
ible objects in the event, by which we mean that the sum on 
individual object-level px values does not include the missing 
transverse energy. The supersymmetric data sets in Figure [1] 
are the total set of all production channels ("All SUSY") as 
well as the subset of these which represent {Ni/Ci) + {q/g) 
associated production ("Semi-Strong"). For the sake of com- 
parison we have generated 5 fb~^ each of tit and QCD di- 
jet production at ^ = 14TeV using PYTHIA 6.4 with 
level one triggers. Distributions for these Standard Model 
background components are also displayed in Figure [T] All 
events have been passed to PGS4 to simulate the detector re- 
sponse. Therefore all objects appearing in Figure [T]should be 
understood as detector objects and not parton-level objects. 
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FIG. 1. Distribution of transverse shape variables tt, and q'rp for 
the bino-like case of point SPS 1 A. In each figure the solid lines rep- 
resent the total SUSY production (thin line) as well as the subset of 
semi-strong SUSY production channels (heavy line). Distributions 
are also given for QCD dijet production and t,f pair production. Each 
plot has been normalized to a constant total number of events. 



To malce comparison of the relative shapes easier the distribu- 
tions have been normalized to constant numbers of events. 

Each variable is effective at separating semi-strong produc- 
tion channels from the overall supersymmetry signal, though 
some of the objects defined in equations (|6]l through ( fTOl i are 
correlated with one another In particular the variables and 
rip are extremely well correlated, though all other pairs of ob- 



servables displayed in Figure[T]are only moderately correlated 
with one another Note that rT,rip — > 1 for a single object 
recoiling against missing energy and vanish for a perfectly 
spherical/isotropic event, or for two antipodal objects in the 
event. 

By contrast, qx 1/2 in the antipodal Umit, 2/3 in the 
isotropic limit, and vanishes in the single object limit. In fact, 
for an N-object, ZAr-symmetric event 



(11) 



QCD di-jet events typically take on values very close to qt = 
1/2, but can often be slightly larger if more than two jets are 
reconstructed in the event. Likewise, distributions for pair 
produced objects, Uke gluinos or tops, tend to center around 
qx = 1/2, but can vary depending on the spread of the event. 
For these reasons, the qx variable is an excellent discriminator 
between semi-strong SUSY events and the backgrounds. 

Note that the variable is the only one of the variables 
in (|6| through (fTOl i defined to explicitly involve the missing 
transverse energy. If the energy and momenta of all visible 
decay products were properly measured we would expect the 
maximum value of to be less than unity. However, the de- 
tector simulator P GS 4 imposes minimum pT requirements for 
object reconstruction and includes the effect of mismeasure- 
ment of object-level pT- Therefore the extracted value of 
can exceed the sum of the transverse momenta of the 'visible' 
objects in the event, and one finds > 1. This is, in fact, 
what happens in the semi-strong processes which interest us; 
we shall return to this issue in the next section. 

The distributions for these shape variables are qualitatively 
very similar for the extreme wino-like scenario. In Figure 
we compare our preferred variables and qx for the bino 
and wino extremes of SPS la. These correspond to the first 
two model lines in Table [V] The agreement is excellent be- 
tween the two cases. Once again all distributions have been 
normalized to equal numbers of events. 



IV. APPLICATIONS OF EVENT SHAPE VARIABLES 

The examples from the previous section clearly indicate 
that the shape variables from Section|lII]are effective at select- 
ing the sub-dominant contributions from semi-strong associ- 
ated production of squarks/gluinos and the lightest neutralino. 
In this section we will outline a procedure for utilizing these 
variables that is sensitive to the wavefunction composition of 
the LSP, gradually adding additional realism as we proceed. 

As demonstrated in the previous section, the variable rip 
defined in (O is one of the more effective shape variables 
in singling-out the semi-strong production processes. It will 
therefore play a central role in the analysis outlined below. In 
Section we indicated that some of these sub-processes will 
have an event shape topology very similar to that of the as- 
sociated production of the LSP neutralino - particularly for 
the wino-like and Higgsino-like extremes. In Figure[3]we plot 
the distribution for the variable rip for the combined semi- 
strong sub-processes and contrast it with the pair production 
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FIG. 2. Distribution of transverse shape variables rr, r'j, and q'rp 
for tlie wino-like case of point SPS lA. Plot labeling is identical to 
Figure [T] 



of strongly-coupled superpartners (gg, gq and qq) and with 
the distribution of rip for the Standard Model background. 
Here we have not taken a sophisticated approach to the back- 
ground generation, but have instead generated 5 fb^^ each of 
t/i and 6/6 pair production, high-p^ QCD dijet production, 
single and Z-boson production, pair production of elec- 
troweak gauge bosons W~, Z and Z Z), and Drell- 
Yan processes. Events were generated at ^/s = 14 TeV using 
PYTHIA 6.4 with level one triggers. 

The region of interest for our purposes is the region where 
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FIG. 3. Distribution of shape variable for SPS lA-based wave- 
function extremes. The bino, wino and Higgsino extremes are given 
in the top, middle and bottom panels, respectively. All data sets in- 
volve a single overall cut of > 175 GeV and the subsequent data 
sets have been normalized to constant numbers of events. 



rip > 1. We immediately observe that the Standard Model dis- 
tribution is relatively flat in rip, but a sizeable fraction of the 
background populates the rip > 1 bins. Clearly, any proce- 
dure for utilizing shape variables to study the wave-function 
of the LSP will require excellent background rejection. We 
will return to this issue shortly. Provided backgrounds can be 
brought under control, it is clear that the majority of the su- 
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persymmetric particle production (represented by the strong- 
strong data sets) can be well distinguished from the semi- 
strong processes. In addition, the shape and peak location of 
these semi-strong distributions (see Figure |2) is clearly influ- 
enced by the wavefunction composition of the LSP. As antic- 
ipated, the strongest signal for semi-strong production in the 
~ 1 region occurs for the wino-like extreme, while the 
weakest signal is the Higgsino-like extreme. For SPS lA we 
expect that examination of events with > 1 will clearly 
distinguish bino-like and wino-like cases from the Higgsino- 
like case. 

The features in these distributions which distinguish be- 
tween the three wavefunction extremes at large values of 
persist even after all supersymmetric production processes are 
included in a single dataset. Figure |4] shows the distribution in 
r'rp for all SUSY processes for each of the benchmark models 
in Table rv] No event selection cuts have been applied to these 
distributions, apart from the level one trigger requirements of 
PGS4, and the distributions are again normalized to constant 
numbers of events. The wino-like extreme always gives the 
largest event rate when the curve is integrated from rip ^ 0.8 
for all four model variants, while the relative sizes of the event 
rates for the bino-like and Higgsino-like extremes will depend 
on the mass ordering between the gluinos and the squarks. We 
note in particular the clear peak in the wino-like distribution 
for all four cases near ~ 1. 

At first glance the results of Figure |4] seem incongruous If 
the property ^ 1 is truly the hallmark of semi-strong pro- 
duction processes, then the simplified discussion in SectionHJ 
would suggest that the bino-like extreme should give a large 
number of events in this region for all benchmark models, 
while the Higgsino-like extreme should be nearly vanishing. 
Some understanding of the puzzle can be obtained by consid- 
ering the relative cross-sections for the different components 
of the total SUSY production for these benchmarks. In Ta- 
ble|VT]we compare the total SUSY production cross-sections 
for the three wave-function extremes of points SPS lA and 
SPS 2. We further subdivide the total into strong-strong 
production, semi-strong production (including all charginos 
and neutralinos) and gaugino-gaugino production of the form 
NiNj, NiCj and CiCj. The strong-strong production cross- 
sections are unaffected by the composition of the neutralinos, 
as one would expect. But this component accounts for the 
bulk of the SUSY production processes only for the bino-like 
extreme for both benchmark points. For the wino-like and 
Higgsino-like cases the pair-production of electroweak gau- 
ginos is a significant component of the total cross-section for 
SPS lA events and is the dominant component of the SPS 2 
total cross-section. 

Clearly, then, the shape variable rip is sensitive not only to 
associated production of strongly-coupled superpartners with 
electroweak gauginos, but also to electroweak gaugino pair 
production events. The 'contamination' is not as severe as 
the numbers in Table IVTl would indicate, however The bulk 
of the gaugino pair production events for the wino-like and 
Higgsino-like extremes involve one or more 'effective LSP' 
states from Table HIH The subsequent decays involve soft lep- 
tons or jets. Therefore many fail to pass the level one trigger 



SPS1A: All SUSY Processes 



— Bino LSP 

— Wino LSP 

— Higgsino LSP 




« 0.035 
<u 

£ 0.030 



SPS1A' : All SUSY Processes 



r:^ (Visible) 



— Bino LSP 

— Wino LSP 

— Higgsino LSP 




SPS2: All SUSY Processes 



r:^ (Visible) 



- Bino LSP 
Wino LSP 

- Higgsino LSP 




1 1.2 

r:^ (Visible) 



SPS2 ': All SUSY Processes 



— Bino LSP 

— Wino LSP 

— Higgsino LSP 




FIG. 4. Distribution of shape variable 
of Table M 



r:^ (Visible) 

for SPS-based benchmarks 



8 



SPS1A: All SUSY Processes (Post Cuts) 
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FIG. 5. Distribution of r'j- shape variable for SPS-based endpoint 
files, after the minimal cuts described in the text. Of particular im- 
portance is the requirement qr < 0.35 which isolates the signal re- 
gion to r'rp ft 0.5. These distributions have been normalized to a 
constant integrated luminosity of 10 fb^'^. 



Extreme 


All SUSY 


SPS lA 
SS SG GG 


All SUSY 


SPS 2 
SS SG GG 


Bino 
Wino 
Hino 


41.48 
63.07 
49.98 


39.99 1.42 0.15 
39.97 2.32 20.88 
39.28 0.33 10.39 


I. 76 
23.87 

II. 29 


1.38 0.04 0.32 
1.38 0.06 22.36 
1.38 0.01 9.85 



TABLE VI. Production cross-sections for benchmark points SPS 1 A 
and SPS 2. In addition to the total supersymmetric production cross- 
section we give breakdowns for strong-strong (SS) production, semi- 
strong production of an electroweak gaugino with a squark or gluino 
(SG) and electroweak gaugino pair production (GG). Cross-sections 
are in units of picobarns. 



requirements. Trigger efficiencies surpass 90% for all strong- 
strong and semi-strong processes, but fall to 7-8% for elec- 
troweak gaugino production for the wino-like and Higgsino- 
like extremes of both SPS lA and SPS 2. After applying 
the trigger efficiencies the electroweak gaugino pair produc- 
tion is still sufficient to upset the simple-minded arguments 
of Section ini but it will not prove fatal to our ability to make 
distinctions between wavefunction extremes using the shape 
variables of Section HIH 

We can use the variable qt, together with kinematic cuts 
on missing transverse energy and transverse mass, to simulta- 
neously diminish the Standard Model backgrounds while iso- 
lating as much as possible the semi-strong component of the 
rip distribution. In our analysis we will impose the follow- 
ing 'minimal' cuts; (1) require > 175 GeV, (2) require 
zero or one isolated lepton, (3) if one isolated lepton, form the 
transverse mass Mt of the lepton with the missing transverse 
energy and require Mt > 125 GcV, (4) require qt < 0.35. 
We plot the distribution in the remaining variable r'rp, after 
these minimal cuts are applied, for our benchmark models in 
Figure |5] These distributions are now normalized to a con- 
stant integrated luminosity of 10 fb^^. The signal region is 
now well-defined as the region with rip > 0.5 and the differ- 
ence between the wave-function extremes is clear for all four 
benchmarks. 

Inclusion of the Standard Model backgrounds, normalized 
to the same integrated luminosity, distorts these distributions, 
even after the imposition of the minimal cuts of the previous 
paragraph. For example, the upper panel of Figure 0adds the 
Standard Model background sample to the SPS lA distribu- 
tions in the top panel of Figure |5] Much of the information 
contained in the shapes near ~ 1 is washed out by the 
large Standard Model contribution. Yet the minimal cuts can 
be easily augmented to recover discrimination in this variable. 
For example, the lower panel of Figure|6]increases the missing 
transverse energy cut to > 250 GcV and require the px of 
the leading jet to satisfy p!^^ > 150 GeV. The excess over the 
Standard Model in the region rip > 0.8 is significant and the 
differences between the wave-function extremes is evident. 

A quantitative measure of the power of the shape-variable 
method to distinguish between various LSP wavefunctions 
can be obtained from integrating the distributions in Figure |6] 
from some minimum value in rip . In Tables IVIII and IVIIII 
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FIG. 6. Distribution of shiape variable for SPS lA with Stan- 
dard Model contribution. The top panel is after the minimal cuts 
described in the text. The shaded region represents the SM back- 
ground; the three curves represent the supersymmetric signal added 
to this background. The bottom panel is after 'augmented' cuts: 
pl^'^ > 150 GeV and fJr > 250 GeV. All distributions are nor- 
malized to a constant integrated luminosity of 10 fb^^. 



r!j. Cut 


Bino 
Events S/VB 


Wino 
Events S/^/B 


Higgsino 
Events S/^/B 


SM 
Events 


r'j. > 0.8 
t-'t > 1.0 
r'r > 1-1 


20,216 61.9 
2,726 14.3 
293 3.5 


24,730 75.8 
5,215 27.3 
667 7.9 


8,229 25.2 
855 4.5 
145 1.7 


150,524 
106,520 
16,760 



r'rp > (r^)inin is clearly capable of distinguishing between 
the three LSP wavefunction extremes in the presence of the 



r'rp Cut 


Bino 
Events S/VB 


Wino 
Events S/VB 


Higgsino 
Events S/VB 


SM 
Events 


r'r > 0.8 
r'j, > 1.0 
r'r > 1-1 


1,083 1.5 
183 0.4 
46 0.2 


15,203 20.8 
5,716 13.4 
979 5.2 


4,949 6.8 
1,371 3.2 
203 1.1 


752,620 
532,600 
83,800 



TABLE VII. Integrated counts and signficances for SPS 1 A. Normal- 
ized to lOfb^V 



TABLE VIII. Integrated counts and signficances for SPS 2. Normal- 
ized to 50 fb"^ 



SM backgrounds. When cuts are optimized it is very likely 
that the goals of our analysis can be achieved for the full LSP 
wavefunction. 



V. CONCLUSION 

Should supersymmetry be discovered in the near future at 
the LHC the energies of the theoretical community will be di- 
rected towards an understanding of the properties of the super- 
partners and the parameters of the underlying supersymmetric 
Lagrangian. Few such quantities are of more general import 
than the wavefunction components of the lightest supersym- 
metric particle. In this paper we have taken a first step towards 
measuring this crucial property using inclusive data collected 
at the LHC. We have introduced new event-shape variables 
and demonstrated their ability to track changes in the contribu- 
tion of sub-dominant SUSY production processes which are, 
in turn, sensitive to the wavefunction of the LSP. This work 
demonstrates the potential power of the technique, which is 
plainly evident from the analysis performed on four simple 
supersymmetric benchmark models. Here we have not tried 
to optimize the cuts imposed, nor sought full experimental re- 
alism. Instead we leave development of a full-fledged analysis 
algorithm to a future work. In particular, it would be of great 
value to analyze how such techniques fare in cases which in- 
terpolate between the wavefunction extremes considered here. 
Ultimately one might anticipate an algorithm that would serve 
as a quantitative measurement of the eigen-components iVi ; of 
the neutralino mixing matrix robust enough to handle a variety 
of supersymmetric models. 



we integrate the tails of the r'rp distribution for three differ- 
ent minimum r'rp values for benchmark SPS lA and SPS 2, 
respectively. For these tables the events were selected after 
applying the minimal cuts described above. For a given over- 
all supersymmetric cross-section - as measured, for example, 
by the inclusive counts of events with fix > 500 GeV or 
some other inclusive variable - the number of events with 
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